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SUMMARY 


The theory of cascades, made -up of a series of flat plates placed 
one behind the other, is extended to the case where the impinging stream 
is not uniform, and the deduced properties of this cascade -flow are then 
applied to the study of the wall interference between such a cascade- 
like bo undary and a vortex-source type of singularity. It is shown that 
the induced velocities, produced by the presence of such a wall, are 
equal to what is obtained by action of a suitably chosen "reflected" 
singularity situated on one side of the wall, together with the action 
of another suitably chosen "transmitted" singularity placed on the other 
side. 


The concepts of a reflection factor and a transmission factor are 
introduced to characterize various types of tunnel boundary, whether 
this is composed of solid wall, open-and-closed sequences, or just a 
free fluid surface . These ideas are then extended to cover the situa- 
tion in which a pair of such walls are allowed to coalesce, expecially 
in the event that one of the walls consists of cpen-and-closed portions 
and the other is a fluid surface._ This latter particular combination 
of boundaries is called a perforated' wall. 

Finally, the interference arising from an arbitrary general singu- 
larity placed symmetrically in between two such perforated walls is 
analyzed, and these results are then applied to the determination of 
the particular kinds of wall geometry that will produce no interference 
effects in the case of a slender lifting wing and. also in the case of 
a symmetric profile having a finite thickness. 


"Studio dell ' Interf erenza delle Gallerie Ae'rodinamiche con Pareti 
a Fessure." Published in Atti della Accademia delle Scienze di Torino, 
vol. 87 , 1952-1953. 
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CHAPTER I 


1. Let attention be directed to a straight wall made up of a series 
of like -sized flat plates, each of chordlength L, and following one 
right after another, with a distance between centers of D. (See 
fig. 1.) 

Let y be the perpendicular distance of any arbitrary general 
point from the straight wall. Furthermore, take Y = y/D as the ratio 
of the normal distance to the "pitch" (or interval between the succes- 
sive repetitions of the gapped wall) . 


The behavior of such a cascade of plates at an infinite distance 
above and below the wall ( i . e . , at the locations where Y = ±<») , that 
is under action of an impinging stream, with a uniform undisturbed 
velocity, denoted by the vector ^ioo, is described by means of the two 


equations 



( 1 . 1 ) 





= V i~ 


+ v 


S°0 


where the velocity below the wall is denoted by the m subscript and the 
velocity in the upper region by the v subscript. Of course, ^Soo 
stands for the induced velocity, produced by action of the gapped wall, 
in a direction parallel with the length of the wall. The magnitude of 
this induced velocity is provided by cascade theory and has the value 


v S co = qy noo (1.2) 

where V noo is the component of which is normal to the plates, 

while q is related to the L/D ratio through the relationship 


q = tan 



(1.3) 


Provided it is understood that the qV noo vector is to be taken 
as lying in the direction parallel with the horizontal extent of the 
cascade-like wall, one may thus rewrite equations (l.l) as 


-¥ 

V m oo 



^ - 
= Via, ” *3Yj3po 

= V leo + qV noo 


(1.4) 
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The sole condition which it is necessary to stipulate in order that 
the equations (1.4) should hold true is that Y = y/D must he infinitely 
large. It may he observed that these equations also remain unaltered 
provided that the L/D ratio stays constant. In view of these facts, 
one may proceed to treat the case where y is merely of finite size, 
however, hy simply letting L and D tend toward zero in magnitude, 
hut in such a manner that the ratio between them is retained at a fixed 
constant value. The assumption will henceforth he explicitly adopted 
that the pitch and the chordlength are, in fact, now both to he infini- 
tesimally small, in comparison with any other of the distances that 
might enter into the argument. In consequence of this assumption, there- 
fore, it follows that equations (1.4) become converted into just 


V m =?i - qV n 

■4 — ) —4 

V v - V ± + q.V n 


(1-5) 


and the gapped wall should thus he looked upon as merely a plate which 
is refracting the stre amlin es . 


On this basis it may also he assumed that the relationships given 
as equations ( 1 . 5 ) will not only hold true for a uniform stream hut that 
they will, in addition, also- he valid for any velocity distribution 
whatsoever, provided merely that the pitch and chordlength are taken 
to remain infinitesimally small with respect to the radii of curvature 
of the streamlines existing in the flow just below and just above the 
perforated wall. 


2 . Now let the flow field he referenced to a complex planar coordinate 
system such that the (reals) x-axis coincides with the direction of the 
straight wall, while the (imaginaries) iy-axis is in the direction lying 
perpendicular to the wall; the origin is taken to he at any arbitrary 
location. (See fig. 1.) The "outside" region of the flow is thus the 
infinite half -plane for which y > 0, while the "inside" portion of the 
flow is the infinite half -plane corresponding to the points for which 
y < 0. Then upon rewriting the vector operations given as equations ( 1 . 5 ) 
as separate scalar equations, it is seen that 


V - V ix - <l v iy and V vx = v ix + 4 v iy] 



( 2 . 1 ) 


and from comparison of the forms of these relationships it follows that 
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V - V ■ ^ v y “ < 2 -5) 

The behavior of the flow in passing through the cascade -like wall 
is thus completely defined by means of equations (2.2) and (2.3). 

The fact that the presence of the cascaded wall does not change 
the normal component of velocity at all is rather surprising. On the 
other hand, it is apparent that the tangential component of the flow 
undergoes a sudden jump as the flow passes the grid pattern of the wall. 
The magnitude of this jump in the tangential velocity component is pro- 
portional to the normal component of the impinging stream. Now this 
change in the tangential velocity-component, brought about by the gapped 
wall, may be interpreted as though it were actually due to the action 
of a distribution of vortices placed along the x-axis and having a 
strength 7, whose magnitude is given by 

7 = v mx ’ Vyx = -2lV y (2.4) 

For the sake of convenience, let a running coordinate along the 
Xrj-axis be denoted by the symbol g. Then the expression for the com- 
plex potential governing the distribution of velocities generated by 
the presence of the gapped wall is readily seen to be 

w = f 2qV y( ! ) log (z - g) dg (2.5) 

2iti J _oo 

It should be pointed out that everything that has been said so far 
is valid only if the flow field is isoenergetic, and it will be taken 
for granted in what follows that this is_ the situation under examina- 
tion, so long as the case of the mixed boundary condition (the perforated 
tunnel) is not being treated. 

In what follows, it will be very convenient to make use of the 
expression for the complex velocity, which may be written down at once 
as 


dw q 

dz iti 



( 2 . 6 ) 


All possible aerodynamic configurations which are met in practice 
can be built up out of the fundamental situation wherein there exists 
in the field of flow one single arbitrary general singularity of the 
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form 


A = Q + ir (2.T) 

(i.e., a combination of a source of strength Q with a vortex of inten- 
sity r) . In what follows it will he sufficient, therefore, to concen- 
trate attention on the study of the flow field resulting from the pres- 
ence of such a singularity. A, placed at the position denoted by -iY 0 
(see fig. l) from the cascade -like wall, and immersed in a stream flowing 
parallel with the x-axis, the magnitude of whose velocity vector is taken 
to be Vjo at infinity. 


The complex velocity corresponding to such a flow field is seen to 
be (with the aid of eq. (2.6)): 


®LL = v i A 1 

dz 00 2at z + iY Q 
in which it is plainly seen that 



Vy(|) 



( 2 . 8 ) 


(2.9) 


holds true. 

3 . Consequently, it is clear that equation (2.8) turns out to be just 
an integral equation in the complex velocity. The solution of this equa- 
tion will be attacked, however, by making -use of an indirect route, 
wherein one relies on the analogy which exists between this case and 
the similar problem with which one is confronted when dealing with a 
solid wall for the boundary (thus the solid-wall example represents the 
limiting case of cascade-type boundary, for which the chordlength of 
each portion of the wall becomes equal to the "pitch", that is to say, 
in the limiting case, one has L/D = 1, and q = co) . 

In this analogous solid-boundary case, it is well known that the 
requisite distribution of vortices lying along the x-axis that gives 
the proper solution to the problem inside_the tunnel (y < 0) is equi- 
valent in all respects to a singularity A (conjugate of A) placed 
at the point +iY 0 (which is the image point of -iY 0 with respect 
to the straight wall) , while the requisite distribution for those points 
lying outside (y > 0) is equivalent to a singularity -A placed at the 
point -IY 0 (which thus cancels the effect of the original singularity 

located there) . 
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Now "this equivalence, which is evinced in "this special case where 
one is dealing with a solid boundary, suggests the possibility immedi- 
ately of extracting therefrom a generalization which will be applicable 
even in the case where one is concerned with the gapped-wall type of 
boundary. In fact the natural thing to do is to take the expression 
for the complex velocity corresponding to a field point lying "inside" 
in the form _ 


=Y m + a. 1 + 41 1 

dz 2it z + iY 0 2it z - iY Q 


( 3 - 1 ) 


wherein A' stands for a suitable singularity located at the point +iY 0 
(the image-point of -iY 0 ) . The value of this latter singularity may be 
determined by invoking the condition that equations (2.8) and (3.1) should 
be identical; that is to say, it must be true that 


dV 

dz 


ay 

dz 


= + 


A_ 1 
2 it z + iY 0 



( 3 - 2 ) 


But it can be seen by reference to equation (3.1) that the value 
of V y (|) is linked to the singularities from equation (2.9) by the 
relation: 


k dz /z-e " \ dz /z=i_ 

1 + L A_ _ — A 1 1 ( V 

ft - 1Y 0 i + iY 0 | + iY 0 | - iY 0 J 

This value of Vy(|) is now to be inserted into the integr al expres- 
sion arising in equation (3.2). The indicated integration is act ually 
carried out merely by integrating along the whole x-axis, entirely within 
the domain of reals. In order to make use of the powerful methods of the 
complex variable, it should be observed, however, that the integrand is 
always going to be an infinitesimal of the second order along the whole 
circle at infinity in the complex plane. Then, if it is formally agreed 
to perform the integration in more generality around a closed contour in 
the complex plane, consisting of the x-axis and the semicircle at infinity 
for values of y > 0, and if the direction of the path of the integration 
is taken to be that of the positive sense for the variable of integra- 
tion, | (which is now considered to be complex) , it follows t hat the 
sought value of the real integral is equivalent to an integration in the 


Vi) - s 


4 rti 




MCA TM 1429 


7 


complex plane, which can be carried out by having recourse to the theo- 
rem of residues, to give 



A A* A 

| - iY 0 | + iY 0 “ | + iY 0 



1 A 1 A' 

2 z - iY 0 " 2 z - iY 0 


( 3 . 4 ) 


Thus, the complex velocity of the flow now under 
written down at once, by referring to equation (5.2), 


examination may be 
as 


= V + ± i + — 3= (A -A') 

dz “ 2 k z + iY 0 2jti z - iY 0 


dz 


(3-5) 


Hence, because of the equality which must exist between the complex 
velocities for the two flows, one is led to the conclusion that 

A> = S (A' - A) = A — 1— (5.6) 

i q - i 

It is obvious,, therefore, that the reflected singularity A', when 
having the value Just determined, will satisfy the requirements imposed 
on the behavior of the flow in consequence of the presence of the gapped 
wall. Hence, the problem which was set has been formally solved in 
closed form. For the sake of convenience, in what follows, the above 
result will be written more compactly as 

A’ = rA (5.7) 

where r is defined as the reflection factor of the wall. Thus it is 
seen that 


r = . jl- aL _ i a_ 

q - i q 2 + 1 1 q 2 + 1 


(3-8) 


This result may be interpreted by means of the statement: The 

"reflected" singularity for a gapped wall is equivalent to the reflec- 
tion factor (complex) multiplied by the complex conjugate of the origi- 
nal singularity. 


If one now goes on to investigate in like manner what happens when 
the field point lies outside of the wall, where y > 0, it is necessary 
to carry out the Integration indicated in equation (5-4), in the complex 
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plane, by going around the closed curve composed of the entire x-axis 
and the semi-circle at infinity lying in the region where y < 0- The 
sense of this integration is opposite to that of | in this instance, 
and the z belongs now to that half of the plane where y > 0. Thus 
the value of the integral becomes now 


r vd) a* a 1 cF 

A ( A* A A' 

_AI_ 

L * - s e f 

j - iY 0 1 + iY 0 I + iY 0 g - iYqJ 

z - | 

. 1 A 

i a' ... i a f 1 q \ 

( 3 - 9 ) 

2 z + iY 0 

2 z + iY 0 2 z + iY 0 \ q + ij 


The complex velocity pertaining to the downstream flow is thus 
obtained by inserting the above expression into equation (2.8), with 
the result that 


dW v ' 

dz 


- v + — - 

“ 2* z + iY 0 


A. __A_ 1 (i . -JL-) 
2it z + iY 0 i \ q. + ±J 


= V„ + — 


2jt z + ±Y r 



-M 

q + ±) 


Voo + 


A i 

2k z + iY 0 


(i 


-r) 


( 3 - 10 ) 


That is to say, viewed from a point outside the gapped wall, the 
field of flow can be looked upon as though it were generated by a single 
singularity A(l - r) located at the point -iY 0 . In analogy to the 
practice common in optics, this will be called a "transmitted" singu- 
larity. The magnitude of this singularity may be interpreted as being 
the product of the magnitude of the original singularity and the trans- 
mission factor t = 1 - F. It then may be noted that 


T « 1 - r = 1 3— = —A— t i — % — - 

9+i q 2 + 1 d q 2 + 1 


(3.H) 


It is worth especial notice to point out the extraordinarily remark 
able fact that each reflection leads once again to -use of the complex 
conjugate of the original singularity, but nothing of a similar nature 
arises when one works with the transmitted singularity. 

Making use of this observation, one may now, easily see the inter- 
connections in the singularities which one obtains when the original 
singularity lies either above or below the gapped wall. To do this, 
let a sketch be prepared that shows the location of the . singularities 
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and the trace of the gapped wall, lying parallel with the x-axis. (See 
fig. 2.) Now Imagine that this scheme of singularity -placement is 
reflected in a line lying parallel with the x-axis. This reflection 
is made in such a way that each quantity will have a complex conjugate, 
as indicated. 

For simplicity's sake the singularities in the flow which have thus 
far been under examination are written with a bar over them. This is 
done in order to obtain a result which will be directly valid for the 
reflected image of the field of flow just about to be studied, wherein, 
for sake of preciseness, it will be assumed that the original singular - 
ity is located at a point lying above the gapped wall. If the original 
singularity lies beneath the gapped wall, it can be seen, by reference 
to the relationships already sketched in figure 2, that the only thing 
which needs to be modified is that the reflection and tranmission fac- 
tors should be taken to be the complex conjugates of the corresponding 
reflection and transmission factors applying in the case where the 
gapped wall is located below the original singularity. 

4. Everything that has so far been derived may be summarized quite 
succintly merely by the statement: The presence of a gapped wall in 

the field of flow coming from the singularity A is equivalent to: 
if the field point lies inside, a "reflected" singularity rA located 
at the image point, while if the field point lies outside, then a flow 
is generated which is equivalent to the presence of a "transmitted" 
singularity tA located at the place where the original singularity 
was situated. 

The factors of reflection and of transmission can thus be written 
as (provided one makes use of a double-level sign, of which the upper 
is to be selected if the singularity lies above the gapped wall, and 
the lower of which applies if the singularity lies below the gapped 
wall) : 


r = 


q. = q. ± i q. 

q. ± 1 q. 2 + 1 ■*- 1 + q. 2 


T = 1 - 


+ i 


q.7i q 2 + 1 i l+q 2 
These factors are linked to each other by the relations 

t = 1 - r 
r = 1 - t 


(4.1) 

(4.2) 

(4.3) 
(*.4) 
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Of course, as was said earlier in Article 1, the value of q is obtained 
from 

1 = tan (| ^ (4.5) 

where L stands for the solid chord length of each open-and-closed piece 
of the gapped wall, and where D represents .the "pitch", or interval 
between sequential pieces of the breached wall. 

Since all possible singularities, either of a concentrated or dis- 
tributed type, can be considered as acting like sources and vortices, 
the results just obtained can be made applicable to all aerodynamic 
singularities . 

As a sort of check, it can be noted, through recourse to equa- 
tion (4.1), that when L/D approaches unity (the solid-wall case) then 
q becomes infinite and r becomes unity too; that is tojsay, when 
L/D = 1, the "reflected" singularity is simply equal to A, regardless 
of what the relative location of the singularity is with respect to the 
wall. On the other hand, the "transmitted" singularity in thiB case 
is tA, where t = 0, and thus this result jibes with the situation 
with which one already is familiar. 

If the L/D ratio tends toward zero, however, which means that 
the wall effectively disappears (but it does not become a fluid sur- 
face, because it is being assumed that the flow is isoenergetic) , then 
q becomes zero, and consequently r = 0 and t = 1. In this case, 
therefore, there is no reflection taking place and the "transmitted" 
singularity is tA = A, and it thus coincides with the original one. 

This is quite obviously the correct result, because if nothing Is pres- 
ent to disturb the flow, it remains the same as it started. Hence, the 
above derived theory appears to be fully justified, and accords with 
known results in the two limiting cases thus tested. 

5. It was explicitly made clear in the above treatment that only the 
behavior of those Interactions was being examined which apply strictly 
to the case of a gapped wall placed in isoenergetic flow. In actuality, 
however, a perforated wind tunnel would require the study of nonisoener- 
getic fields. This is so because the real perforated tunnel can be 
thought of as being a superposition of a. free fluid surface (that Is, 
this boundary is just the line of demarcation between the tunnel jet 
and the external atmosphere) upon the gapped-wall type of flow. When 
two such boundaries are allowed to coalesce, they constitute together 
what will now be called a perforated -wall tunnel, or simply a perforated 
wind-tunnel. 
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It is well known that the effect of a free fluid surface (two- 
dimensional) operating in conjunction with a singularity A is equiva- 
lent to the presence of a singularity -A placed at the image point 
of A; that is, the reflection factor for a free fluid surface is 

r' = -1 (5.1) 


and consequently 


t’ = 1 - r’ = 2 (5.2) 

In order to study most conveniently the case of such mixed "bound- 
aries, it is "best at first merely to handle each wall separately and 
then let them coalesce. By aid of this process of "bringing the two 
walls into juxtaposition, one may thus deduce what the appropriate 
reflection factor should be for such "mixed" walls. 

To proceed with this analysis of the double walled tunnel, it is 
convenient to denote the gapped wall as P B and the free fluid surface 
as Pf. Also let A = Q + r/i represent a singularity located inside 
the wall. (See fig. 5.) 

It should be recalled that the selfsame wall will have reflection 
and transmission factors of r and t, respectively, or of r and t 
depending upon whether the radius vector, which originates at the singu- 
larity which is giving the reaction at the wall, approaches the wall in 
question from above or from below. 

If one imagines that the radius vector emanating from the singu- 
larity A has the 6ame attributes as would appertain to a ray of light 
issuing from this source, it is quite easy to find the locations of the 
multiple reflections R s , Rf^, • • • which are produced by action 

of the two walls under study, and at the same time it will be equally 
easy to establish where the points of incidence are located depending 
upon whether the ray approaches them from below or above. 

By use of a system of notation which is perfectly obvious, it may 
be stated that the first reflection point is thus given by the relation 

Rg = rgA (5*3) 

In order to arrive at the point fj_, the ray u must undergo a 

process consisting of a transmission at P B characterized by the fac- 
tor a reflection at characterized by a reflection factor r^j 

and in turn another transmission at P s characterized by the factor t b . 
As a result, the ray in question now appears to issue from the reflection 
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point 

as 



and the strength of this reflected singularity is thus given 


Rf-L = (Ar s )r f T s = Ax f r a ‘ 


(5^) 


In order to arrive at the point fg, the ray must undergo a process 


consisting of three reflections and two transmissions. As a result the 
ray will then appear to issue from the reflection point Rfg> and the 


strength of this singularity will be found, according to a process which 
is entirely analogous to what was done in the previous case, to be 


R*>_ = 


jAT s )r f 


f T s r f = Ar f T B 2 r f f s = B^r^. 


( 5 - 5 ) 


The simple repetition of this process of tracing out the paths of 
the light ray thus leads one to write in general: 


Rf (m-1) = Rf n rfrs = R fi( r f r s) n " A^VC^s) 


TT _ 2. - s n 


( 5 - 6 ) 


It is now time to see what happens when the two walls and P s 

are allowed to come closer and closer until they finally coincide. It 
so happens that now all the "reflected" singularities once again come 
together and coalesce in the image point of A, while their magnitudes 
are all gathered into one single "reflected" singularity A*, the mag- 
nitude of which is 


A* = A 


00 

I 


r s + > r f T s 2 ( r f r s) n| 

ne=0 


» Ar. 


m 


( 5 - 7 ) 


where r m stands for the reflection factor of the perforated (or mixed) 

wall. By carrying out the indicated sum of the geometric series, it 
follows that (provided it is taken for granted that t s = 1 - r s ) : 


r -r ir, t 1 ~ 

” * f 1 . rf - 


( 5 - 8 ) 


^The bar placed over the symbols in parentheses is used to indicate 
that one should take the complex conjugate of the indicated quantity. 
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In tlie case in which interest is now focussed, the wall P f is to 
he taken as a free fluid surface, and thus the reflection factor which 
pertains to a perforated tunnel is obtained by making use of equation ( 5 . 8 ) 
(which is true in general) by particularizing it to the case wherein r^ 

is set equal to r f , and where r f = - 1 , and where, besides, r B = ■ % — . 
Thus now 


r - (1 - r s f 3r, - 1 21 2 - - 1 

r m - r s r: = : : 

1 + r s r s + 1 2 q 2 - i q + 1 


(5-9) 


It is also clear that if the perforated tunnel wall happens to be 
below the si n gula r ity, the reflection factor is the conjugate of the 
one just now deduced. 

One may check the above-derived result by letting the L/D ratio 
vary from one limit value to the other and noting what transpires. For 
instance, let L/D = 1 first, and then q = oo (a solid wall and a free 
surface coming into coincidence) and thus the reflection factor in this 
case is r m = 1. On the other hand, if L/D is allowed to be zero, 
then q = 0 (the solid wall disappears but a single free fluid surface 
remains), and in this case the reflection factor is r m = -1. These 
results are just what one would expect. 


It is worthwhile pointing out that 
may write 


L m 


= e 


Id 


r m = 1 , and therefore one 


(5-10) 


where 


d = argument (rp) = G(L/D) 


(5-11) 


The value of this functional relationship ■d = g(l/d) has been 
computed for a number of values of L/D, and the results are reported 
in the appended table for the case where the tunnel wall lies under- 
neath the singularity A. 


l/d 

d 

L/D 

d 

0 

-5.14 

0.5 

-O.93 

.05 

- 3.125 

.6 

-.722 

,1 

-2.54 

• 7 

-.500 

.2 

-1.995 

.8 

-.524 

• 3 

-1.55 

• 9 

-.158 

.4 

- 1.21 

1.0 

.00 
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CHAPTER II 


1. In the next step towards understanding the action of a perforated 
wind-tunnel, it is necessary to take into consideration the effect of 
both of the perforated walls P and Q, one of which lies above and the 
other lies below the singularity A. The distance between these walls 
is assumed to be the height h, and they are symmetrically placed with 
respect to the singularity, which is itself considered to lie at the 
origin of the complex coordinate system to which the flow is referred. 

It is clear that the usual multiple reflections also take place 
here, and in consequence it will be true that the velocity at each point 
in the region of the flow comprised between the two walls will be composed 
of the sum of the undisturbed velocity V*, (the velocity vectors produced 

by direct action of the singularity A, together with the velocities 
arising from interference from all the reflected singularities). 

Now let rp and r^ be the reflection factors for the lower wall 
and upper wall, respectively. Then if one examines the ray Up emanating 

from A and progressing downward from A (this path is denoted in fig. 4 
by means of a solid line) it will be seen to be reflected at the wall P, 
and because of this it will appear to have issued from the singularity 
Pq_ located at a vertical distance downward of -ih. The strength of 
this singularity is given by the relation 



( 1 . 1 ) 


Continuing on with the tracing out of the path followed by the 
ray Up, it will be recognized that it is again reflected from wall Q, 
and it will thus appear to have come from the singularity P-^' , located 

at a vertical height of 2ih, and having a strength given by 


V - p i r q 


= Ar. 


p/ q, 


= Ar P r a 


( 1 . 2 ) 


It is intuitively obvious that the path of the ray Up will be 

periodically repeated in such a way that any two successive "reflected" 
singularities P n ' and P^j. will have the same relation between each 

other as the singularities P^ 1 and A have. Thus it will be true in 
general that 


p» 

n+1 


= P 'r r 
n p q. 


(1.3) 
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and consequently it will be true that 

V =A ( ? p r q) n (lA) 

while at the same time the ordinate location of the singularity V 
will be given as 


iy p • = i2nh (1.5) 

r n 

Returning now once again to take up the question of the continuation 
of the path of the ray Up after it has undergone the last reflection 

considered previously, it will be seen that it is once again reflected 
at the wall P. In consequence of this reflection the ray now will 
appear to be emanating from the singularity P£ which lies at a distance 
-i3h below the wall, and the strength of which is given by 


P. = P, *r = 5r_r r_ = P n r_r_ (l.6) 

2 1 P P9PI9P v ' 

It will again be abund an tly clear, therefore, that any two successive 
"reflected" singularities p n+l and P n will be interrelated in exactly 

the same way that the singularities P 2 and Pp were found to be 
connected above. Thus, in this case, it is obvious that 

Pn+l = P n r q r p (l»7) 

and consequently it will be true that 

P n = P l( F q r p) 11-1 * kp( ? q r p) n ” 1 ^ • 8 ) 

and the ordinate distance of P n is given as 

iyp = -i(2n - l)h (1.9) 

r n 


Py means of this procedure half of the reflections of the singularity 
at A have now been determined. In order to find the other half one may 
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proceed to follow an entirely analogous development which stems from the 
tracing out of the path followed by a ray u^. The route taken by this 

ray will be entirely symmetric to the one traversed by Up (this new 

path is also drawn in on fig. 4) except that now the order in which the 
walls are encountered by the ray Uq will be inverted from what was 

true for the Up ray. It is quite evidently sufficient, therefore, 
merely to interchange the symbols rp and r<j in the above -written 
formulas in order to arrive at the correct expressions for the magnitudes 
of the "reflected" singularities ^ and Q^, which lie at locations 

which are symmetrically placed with respect to the analogous singular- 
ities P n and P n ’ . The sought values appertaining to this half of the 

"reflected" singularities are thus obtained as: 



9n - 

(1.10) 


iy 0 = i(2n - l)h 

(1.11) 


<*»’ = -“(Vp ) 11 

(1.12) 

and 

ly<j n ' - -12nh 

(1.15) 


2. It is now convenient to introduce the factor f(z-z Q ), which is the 

* 

amount by which the magnitude of an arbitrary singularity located at the 
point z Q must be multiplied in order to determine the magnitude of the 
complex velocity produced at the point z. ly making use of this conven- 
tion and by having recourse to the expressions just derived in the pre- 
vious paragraph, one may write down the value of the complex velocity as 
a sum of all those complex velocities which are produced by action of the 
multiple reflections in the two perforated walls. This total value for 

the complex velocity is thus 


dw 

dz 


AJ 


00 

X r p(f ( l r p) n ~ 1 f t + ±(2n “ l)h J 


n=l 


oo 

Y r i( ? p r i) n_lf & ■ i ^ 2n ■ i * ^ h ] 


n=l 
A 


? + 


00 oo 

Y ( 5 P r< l) n:f ^ Z ■ i2nh ) + Y ( 5r< 3- r p) nf ^ Z + i2nh ^ 


n=l 


n=l 


( 2 . 1 ) 
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Since in actual practice it will be tame that the most commonly met 
tunnel configuration will be one for which both of the perforated walls 
have the same reflection properties, and since there will be a good deal 
of simplification resulting from such an assumption, it will be assumed 
that this is the type of tunnel with which the following discussion will 
be cone earned. In this case, it is true that 


r p = (2.2) 

and consequently, equation (2.1) will become simplified to Just 


dv 

dz 


00 



+ 


[n=l 


2^ (rp) 2rL “ 1 f jz - i( 2 n - 
n=l 



+ 



r p 2n f(z + i2hh) + 


n=l 


i2nh) > 


( 2 . 3 ) 


One can bring about another great simplification in the work if 
the complex velocity is only sought for points lying along the axis of 
reals (the x-axis) . If this limitation is agreed upon, then it will 
be seen that the arguments of f are complex conjugates in pairs, and 
likewise the sums of such pairs thus also will be complex conjugates. 
Consequently, the Siam total of the arguments is juBt twice the value 
of the real part of each of the individual ones, and it follows that 


00 00 
= 2AR r p 2 * 1 - 1 * [x + i(2n - l)h] + 2AB. ^ r^fU + i2nh) 

' Z yi n=l n=l 


00 

= 2R(A)R 


n=l 


fp 211 ^f (x + i(2n - l)h| + r p 2n f(x + i2nh)| + 


121 ( A) R ^ ■Jr p 2n f(x + i2nh) - rp 211 ' 1 ^ + i(2n - l)hj| (2.4) 
n=l *• 
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It may be noted that each term of the Indicated sums is the sum 
of two terms which involve the quantities 2n - 1 and 2n. It is con- 
venient to reletter these multipliers as m and m + 1 , respectively. 
Hence one may now rewrite the expression for the complex velocity more 
concisely as 


00 OP 

= 2R(A)R *p m f(x + + 2il(A) R Y (- r p) mf ^ x + ^ (2.5) 


m=l 


m=l 


Regardless of what the expression for the function f happens to 
be, it is worthwhile pointing out that it is always possible to develop 
equation ( 2 . 5 ) in a MacLaurin series provided one only moves away from 
the origin a short distance, that is to say, provided the value of x 
is small in comparison with the width of the tunnel, h. Upon carrying 
out such a development one obtains the complex velocity as 


(£) = Yj r p mf ( imh) + 

m=l 


2il(A)R Y 


m=l 


^-r p j m f (imh) + 


X 


00 r~ 

2R(A)R V* rJ 11 — f(x + imh)l 

Z_i P L^x J x=0 


m=l 


2il(A)R £ (- r p)”fe f <* + “ 

m=l „ 


+ . . . 


In order to carry on the analysis any further it is necessary to 
select a specific functional relationship for f(x + imh) for each 
particular case, depending upon what type of singularity is under study. 

5 . In the case where the singularity is made up of a source of strength, 
Q, and a vortex having a circulation, P, that is, for the case where 



f(x + imh) = — - — ( 3 - 1 ) 

2 it(x + imh) 


and consequently 
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m=l m=l 


x 

h 



+ 


(3.2) 


provided the series development is carried out only as far as the first 
term in x. 


4 . In the case where the singularity is a doublet, of moment M, with 
its axis orientated so as to be parallel with the x-axis, the appro- 
priate form for f is 

f(x + imh) = .+ i (4.1) 

2jc(x + imh) 2 


and consequently 



Ot.2) 


provided the series development is once again merely carried out as far 
as the linear term in x. 

5- In the preceding two Articles, one is confronted with the necessity 
of summing certain series which are expressed as follows: 

oo / , \ nL co 

£ V l_gl = a y (±l) m cos ± 1 sin mfl ( 5>1 ) 

4 (im)* 4 (im) t 

m=l ' m=l 

where t = 1, 2, 3, • • • . 

The latter form for the series results from acknowledgment of the 
fact that 


r-. = — - - = cos •fl + i sin a (5-2) 

* 2q - i + 1 

inasmuch as |rp| =1. While it is also true that the restriction 

-n < a = arg fr^\ < 0 


holds. 


( 5 - 3 ) 
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The case a = 0 corresponds to that of a solid wall boundary, while 
the case a = -jr corresponds to the situation where the boundary is a 
free fluid surface. 


Now, depending on whether one employs the + sign or the - sign in 
equation (5.1) , different trigonometric series are obtained, and like- 
wise when different values of t are inserted into the formulas differ- 
ent trigonometric expressions are generated; but these are all easily 
summable, and the several cases of interest are presented below: 


Case where t = 1 and the sign is 


oo -m 

V -2— 

00 

Y 1 sin ma 

_ * 

_ a 

La im 

L m 

I 

2 

m=l 

m=l 




(5.4) 


for a lying within the open interval running from 0 to -2«, while 


*1 


m 


iS-mO 

in 


for a = 0 or a = - 2 k 


Case where t = 2 and the sign is + : 



m=l 


Z cos ma 

m2 



(5.5) 


for a lying anywhere on the closed interval from 0 to - 2 k , including 
the end-points. 


Case where t = 5 and the sign is + : 


pj Y r p m _ Y sin ma _ «a^ + a£ 

L> laP L-i m^ ^ 1-2 


|2 . & , j^a 

~ 


(5-6) 


m=l m=l 

for a lying anywhere in the whole closed interval from 0 to -2k. 
Case where t = 1 and the sign is : 


* I - X 


m 00 


sin ma 


1 

2 


m=l 


m=l 


m 


(5-7) 
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for d lying within the open interval running from -it to +it while 



m=l 


0 for d = ±it 


Case where t = 2 and the sign is 

(5 ' 8) 

m=l m=l 

for d lying anywhere on the closed interval from -it to +it, including 
the end-points. 


These expressions are all that is required in order to be able to 
describe the behavior of the velocities along the axis of symmetry of 
the wind tunnel in the neighborhood of the singularity. By substitution 
of the above -evaluated series in the expression given as equation ( 3 . 2 ) , 
it is seen that the sought formula for the complex velocity along the 
axis of symmetry turns out to be, in the case of a singularity of the 

form A = Q + ^: 

{S') = --S-U + i) -_£_* + 

\dz y x £jtll 2inh 

X 

h 


a 

2 ith 



(5.9) 


while if the properly evaluated trigonometric series are substituted 
into the expression given as equation (4.1) it is seen that the sought 
formula for the complex velocity along the axis of symmetry turns out 
to be, in the case of a doublet with moment M: 



M 

2 

2 itn 



2 x M ( it 2 d ltd 2 d^\ 


( 5 . 10 ) 
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CHAPTER III 


1. After having gathered together all the necessary "basic tools as 
derived above, it is possible now to go ahead with a study as to just 
how one might eliminate the interference effect at certain important 
locations in a perforated wind tunnel. 

In this study two important limiting cases will be examined; viz., 
in one case the infinitesimally thin lifting wing will be treated (the 
antisymmetric problem), while in the other case the symmetric wing of 
finite thickness but at a zero angle of attack, will be studied (the 
symmetric problem) . 

As is’ well known, the lifting wing (antisymmetric case) is handled 
by imagining that a vortex is concentrated at the quarter chord point 
of the airfoil with chordlength l, and the boundary condition that must 
be satisfied is that the streamlines become tangent to the mean camber 
line of the profile at the three-quarter chord point. 

Now, in general, because of the interference effect of the surrounding 
wind tunnel, it will be true that the local angle of attack at the three- 
quarter chord point will no longer be_the same as that angle given by the 
difference in direction between the V^, direction and the chord line of 

the profile, but instead this local angle of attack is increased by the 
amount 


e 



( 1 . 1 ) 


where Vy is the component, taken in the direction of the y-axis, of the 
velocity induced at this point by interference action of the tunnel walls - 

Consequently the boundary condition which now has to be satisfied 
is 

r = Jtz(a + e)V m (1.2) 

where r is the circulation existing around the profile in question, and, 
of course, a is the angle of attack of the chordline of the airfoil with 
respect to the free-stream direction, and I represents the chordlength. 

If one wishes to eliminate the wall interference it is obvious that 
this means that one must counteract or cancel out the increment e, and 
thus the Vy component must be annulled at the three-quarter chord point. 
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Let it now "be assumed that the profile is located in the tunnel in 
such a way that the point where the vortex r is considered to he concen- 
trated lies on the axis of the tunnel (the x-axis) and that this point 
is also the origin of coordinates. Thus the three-quarter chord point 
then lies at a distance x = l/2 downstream 


By reference to equation (5-9) of Chapter II, it will he seen that 
under these stipulations one may write 


/dw\ r 


W^i/a 

2h\2 6/j 


(1.3) 


and thus the vertical component of induced velocity may he eliminated hy 
making 



1 ± 


1 + 


«2 
h2 12 


(l.fc) 


Since 9 = arg(rp), it is less than zero and thus one must select 
the negative sign in the above formula. 


On the other hand, if one would rather consider that 9 is the 
fixed quantity, then it will he necessary to choose the ij h ratio so 
that 


if h 


129 



(1.5) 


2. Now the symmetric problem is to he treated. 

It is well known in this case that the symmetric profile can he 
simulated hy employing a distribution of sources and sinks whose total 
strength is actually zero, and this distribution is equivalent (as far 
as its effects on the field of flow, at points far away from the loca- 
tion of these singularities, are concerned) to the effect of placing a 
doublet, -whose moment is going to he proportional to the cross-sectional 
area of the profile, i.e., M = -V^S, at the centroid of this area. 


Let it he assumed, therefore, that the centroid of the profile is 
placed so as to coincide with the origin of the coordinate system for 
the tunnel flow, and let x he the abscissa val u e for the point at 
which one wants to cancel out any interference effect that would arise 
from action of the constraining wind tunnel walls. The magnitude of 
this induced velocity will thus he expressible as 
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M 

2kh 2 



+ It'S + 


_ 2x / k 2 fl + + 

2) h V 3 2 T/ 


( 2 . 1 ) 


and In order to cancel It, one most obviously impose the condition, there- 
fore, that (since d is to take on only negative values): 


*£ + _ 2x(k£i + %£. + 

3 2 h \ 3 2 6/ 


0 


( 2 . 2 ) 


This relationship may he rearranged into a more 
work with, and upon carrying out this simplification 


convenient form to 
it will he seen that 


x = + 3 (* + 0.422k) (d + 1.578k) 
h 2 a(a + jt)(d + 2k) 


If the value of ■d is confined to the closed interval running 
between zero and -k, then the ratio x/h varies continuously from a 
value of +oo to a value of -oo and passes through the point zero when 
■d = -0.422k. 


The behavior of this relationship for x/h is of particular inter- 
est in the neighborhood of the value where d = -0.422k (or where 
x/h = 0, that is). In the region of the ■fl close to i = -0.422k the 
expression for x/h may be developed in a series to give simply, to 
first order, that 

| = -0.462(« + 1.323) + . . . 

where " (2.4) 

$ < 0 

Conversely, it is easily seen what the expression for d will be 
as a function of x/h in this case, and this relationship may thus be 
written as 

a = -*(2.016) - 1.323 (2.5 

In the case where an airfoil has central symmetry (both midchord 
as well as fore-and-aft symmetry) it is obvious that the interference 
from the wall will turn out to be zero at the location of the center of 
symmetry, provided the value of the porosity is such that d = -1.325 • 


Translated by R. H. Cramer 

Cornell Aeronautical Laboratory, Inc., 

Buffalo, New York 
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(a) Field point inside. 

Figure 2 


mCA TM 1429 


• Z 


• T A 

Axis of 
reflection 

• rA 


• z 

(b) Field point outside. 






